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ABSTRACT 

' We present the initial results from our search for high redshift, z > 6, quasars using 

\ near infrared data from the UKIRT Infrared Deep Sky Survey (UKIDSS) Large Area 

7— I i Survey (LAS). Our analysis of 106 deg 2 of sky from Data Release 1 (DR1) has resulted 

(N ; in the discovery of ULAS J020332.38+001229.2, a luminous (Jab = 20.0, J Vcg a = 

^ — 1 , 19.1, M1450 = —26.2) quasar at z = 5.86. Following candidate selection from the 

combined IR and optical catalogue data and stacking of multiple epoch Sloan Digital 
Sky Survey (SDSS) data, we have obtained optical spectroscopy for the only two 
43 z > 6 quasar candidates. The VLT FORS2 spectrum of ULAS J020332.38+001229.2 

Qh shows broad Lya + Nv 1240 emission at ^8350 A and an abrupt continuum break 

1 , due to absorption by the Lya forest. The quasar is not present in the SDSS DR5 

■ catalogue and the continuum spectral index of a=-1.4 (F w oc v a ) is redder than a 

\ composite of SDSS quasars at similar redshifts (a=-0.5). The discovery of one z ~ 6 

C^J quasar in ^100 deg 2 in a complete sample within our selection criteria down to a 

k* ' median depth of Yab = 20.4 (7c) is consistent with existing SDSS results. We describe 

. £h ! our survey methodology, including the use of optical data from the SDSS and the 

. highly effective procedures developed to isolate the very small surface density of high- 

' probability quasar candidates. 



Key words: cosmology: observations, galaxies: quasars: general, galaxies: quasars: 
individual: ULAS J020332. 38+001229.2 



1 INTRODUCTION 

Quasars are powerful probes of the high-redshift Universe 
since, by virtue of their high luminosities, they can be ob- 
served out to immense distances with large look back times. 
The host galaxies of the highest redshift quasars are proba- 
bly still forming and most likely occur in overdense regions in 
the matter distribution of the early Universe. Studies of their 
space densities, host galaxies and local environments provide 
potentially powerful constraints on theories of the formation 
and growth of supermassive black holes, the growth of large 
scale structure and the formation and evolution of the first 
galaxies (e.g. Kauffmann & Haehnelt 2000). 

In addition to being of intrinsic interest, bright high- 
redshift quasars are particularly valuable as probes of the 
Universe via absorption studies of cosmologically distributed 
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intervening material and can be used to determine the bary- 
onic content and physical conditions (metallicity, tempera- 
ture, ionization state) of the intergalactic medium (IGM) 
and gas within intervening galaxies. 

Finding quasars above redshift z > 6.0 is of crucial 
importance since it is in the epoch of reionization between 
the redshift where the Stromgren spheres of the reionization 
sources overlap (zovl = 6.1 ± 0.15, Gnedin & Fan 2006), 
as shown by the observations of the Gunn-Peterson optical 
depth at z ~ 6, and the lower bound (1<t) to the epoch of 
reionization from the 3-year WMAP polarization results at 
z — 8.2 (Page et al. 2006) which corresponds roughly to the 
epoch when the Universe was half-ionized (by volume). 

Over the last few years there has been considerable 
progress in wide area searches for high-redshift quasars. 
Storrie-Lombardi et al. (2001) reported the final results of 
a large survey for z ~ 4 quasars over ~8000deg 2 which re- 
sulted in ~50 quasars in the redshift range 4.0 < z < 4.8. 
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This survey used photographic data and was based on the 
drop across Lya caused by the Lya forest, depressing the ob- 
served flux in the _B-band relative to the J?-band. The advent 
of the CCD-based Sloan Digital Sky Survey (SDSS, York et 
al. 2000) extended wide-area surveys to longer wavelengths 
with the inclusion of the i- and z-bands, enabling searches 
for quasars based on i — z colour and redshifts of z ~ 6 were 
reached (Fan et al. 2000). There are now around 20 z ~ 6 
quasars known, with the highest redshift quasar at z = 6.42 
(Fan et al. 2003, 2006). 

However, finding quasars beyond z — 6.4 using the 
SDSS is practically almost impossible because the sources 
become too faint in the reddest passband (the z-band) due 
to absorption by the intervening Lya forest. To reach higher 
redshift requires wide-field near-infrared imaging to sub- 
stantially fainter fluxes then the 2 Micron All Sky Survey 
(2MASS, Skrutskie et al. 2006). This is the parameter space 
explored by the UK Infrared Telescope (UKIRT) Infrared 
Deep Sky Survey (UKIDSS) Large Area Survey (LAS). One 
of the main science drivers of the UKIDSS LAS is the dis- 
covery of quasars with z > 6. 

In this Letter, we report the first high-redshift quasar 
discovered with UKIDSS, ULAS J020332.38+001229.2 at 
z = 5.86. All magnitudes are given in the AB system using 
the Vega to AB conversions provided in Hewett et al. (2006) 
and a A-dominated cosmology with Ho = 70 kms -1 Mpc -1 , 
Om = 0.3, and f2 A = 0.7 is adopted. 



2 SURVEY IMAGING DATA 

The near infrared data are from the UKIDSS LAS (Lawrence 
et al. 2006) which is being carried out using the UKIRT 
Wide Field Camera (WFCAM, Casali et al. 2006) on the 
3.8-m UKIRT. WFCAM has a sparse packed mosaic of four 
2048x2048 pixel Rockwell Hawaii-II arrays with a pixel scale 
of 0'.'4 pixel -1 . The instantaneous field of view is 0.21 deg 2 . 
The arrays are spaced by 0.94 times the detector width. 
A mosaic of four pointings gives a tile with a continuous 
covered area of 0.77 deg 2 . See Dye et al. (2006) for technical 
details of UKIDSS. 

The LAS consists of nominal 40-s exposures in the Y-, 
J-, H- and A'-bands. The Y-band (0.97-1.07 /im) is specially 
designed to fall between the z-band and the J-band, occu- 
pying the clean wavelength range between the atmospheric 
absorption bands at 0.95 /im and 1.14 pm. The Y-band is 
needed to distinguish high-redshift quasars from the more 
numerous cool stars, such as L- and T-dwarfs (Hewett et al. 
2006, see also Fig. [[J. 

In this Letter, we use data from the Data Release 1 
(DR1, Warren et al. 2006). The UKIDSS observations were 
obtained during the period 2005 August to 2006 January. 
The data was released to the ESO community in 2006 July 
(see Warren et al. 2006 for further details). The median 
5<t point source depths in AB magnitudes are Y\b=20.8, 
Jab=20.5, #ab=20.2 and K A B=20.1 (Warren et al. 2006). 
We have restricted our analysis to the UKIDSS LAS area 
that is imaged at multiple epochs by the SDSS Southern 
Survey (-25° < RA < 60°; -1.25° < Dec < 1.27°, called 
Stripe 82 in the SDSS footprint, York et al. 2000). This re- 
gion of sky is being imaged by the SDSS collaboration in 
5 passbands (u, g, r, i, z) repeatedly during three months 




Y-J [AB] 

Figure 1. i — Y vs Y — J diagram in the AB magnitude sys- 
tem illustrating simulated colours of stars, elliptical galaxies and 
quasars (from Hewett et al. 2006) and observed colours from 
UKIDSS and SDSS. Also plotted are data from 20 deg 2 from the 
UKIDSS Early Data Release. The cloud of blue points are objects 
classified as stars, red points are objects classified as galaxies. The 
region in colour-colour space occupied by high-redshift quasars is 
outlined by the dashed lines. Quasar ULAS J0203+0012 is plotted 
as open circle with error bar (SDSS DR5 j-band photometry) and 
as solid circle (co-added i-band photometry). The solid square is 
the M dwarf ULAS J022905.12+005506.5. 



(September, October and November) in each of three years 
(2005-2007) as part of the SDSS-II project. In addition 
Stripe 82 was imaged 10-20 times during the SDSS-I project 
(see Adelman- McCarthy et al. 2006). The SDSS data is re- 
leased as DRSN1 (Sako et al. 2005) and DRsup (Adelman- 
McCarthy et al. 2006). The area covered by the LAS in both 
Y- and J-bands in the SDSS Southern Survey area is 106 
deg 2 . In this area ~950 000 objects are detected in both Y 
and J with a signal-to-noise ratio in Y ^ 7. 

We matched the UKIDSS DR1 data with the catalogues 
from the SDSS DR5, using a search radius of 2 arcsec. The 
colours of objects are computed using the SDSS PSF mag- 
nitudes and the UKIDSS apermagS magnitudes which are 
computed using a diameter of 2 arcsec, corrected for aper- 
tures losses on a pointing by pointing basis. 



3 QUASAR SELECTION 

Candidate high-redshift (z > 5.8) quasars down to a signal- 
to-noise ratio of ~7 in Y, corresponding to a median depth 
of Yab = 20.4, were selected on the basis of their blue Y — J 
UKIDSS colour and, either, (i) red i — Y colours from our 
combined UKIDSS-SDSS catalogue, or, (ii) absence in the 
SDSS i- and/or z-band catalogue. The colour-colour dia- 
gram shown in Fig. [T] demonstrates the principles behind 
our selection criteria. The colour cuts used were i — Y > 2.5 
and Y — J < 0.6. The Y — J cut is used to avoid the L and 
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T dwarfs. The i — Y cut selects quasars above z = 5.8 and 
discriminates against foreground Galactic stars. 

An additional requirement for candidate quasars was 
that they should be undetected in the SDSS it-, g- and r- 
bands. We regarded an object as undetected in SDSS, if the 
SDSS PSF-magnitude was fainter than the 3cr limiting mag- 
nitude of the SDSS 2048 x 1489 pixel (13'5 x 9'.8) field in 
which the object was found. We made an empirical determi- 
nation of the 3ct magnitude limit for each SDSS field from 
the SDSS catalogue. If an object was detected in the i-band, 
we also required that i — z > 2.2 following the criterion used 
by Fan et al. (2003). 

We performed various steps to remove false positives 
from our candidate quasar sample. These steps included: 



(i) Removal of intra-quadrant cross-talk images. Cross- 
talk images are artifacts caused by bright stars (see Dye et 
al. 2006 for a detailed description). The cross-talk images 
appear offset from the bright star by multiples of 51.2 arc- 
sec (or 128 pixels with scale tf.'A) corresponding to a single 
read out channel on the same quadrant. More than 99 per 
cent of the z-dropouts are cross-talk images. We used the 
2MASS catalog to locate all stars in the surveyed area down 
to Jab = 15.5. We then used the positions of these stars 
to remove all candidates that were located within a radius 
of 5arcsec of multiples of 51.2 arcsec in the North/South 
or East/ West direction as appropriate for the quadrant (see 
Dye et al. 2006) out to ~6 arcmin on the F-band images 
from the 2MASS stars. 

(ii) Independent repeat photometry on the SDSS images 
using the pixel data of the SDSS to measure the photom- 
etry of candidate quasars in the g-, i- and z-bands. This 
removed objects that do not appear in the SDSS catalogues 
but are present on the images, such as in the case where 
objects lie near perturbing brighter objects and the edges of 
the SDSS data. This photometric analysis also provided 3a 
upper limits. 

(iii) Co-adding of up to seven epochs of SDSS data from 
DRsup and DRSN1 taken in photometric conditions to de- 
tect objects that lie below the SDSS single-epoch catalogue 
limits in the optical i- and z-bands. The co-added images 
were obtained by averaging multiple single-epoch SDSS im- 
ages centered on the location of the UKIDSS object, scaled 
to a common zeropoint using the DR5 catalogue. 

(iv) Screening by eye to identify spurious candidates 
which had not been removed at any earlier stage (e.g. mov- 
ing objects such as asteroids and corrupted objects near per- 
turbing bright objects). UKIDSS observations in each wave- 
band are not simultaneous and the time interval between 
consecutive observations in different wavebands is generally 
in the range 30 to 120 mins. Asteroids could be recognized as 
they move between the observations in the different UKIDSS 
bands. An object that appears on more than one band and 
shifts by > 1 arcsec hr -1 could be identified as moving, down 
to the flux limit we considered. 



These steps reduced the number of candidates from 
>5000, the majority of which were cross-talk artifacts, down 
to 10, the follow-up observations of which are described next. 



Table 1. Photometry of the z = 5.86 quasar, ULAS 
J020332.38+001229.2. Magnitudes are in the AB system. 



Band 


Source 


mag [AB] 


Date 


i 


SDSS DR5 


> 23.1 


2002-09-05 




Co-add (DRSN1, DRsup) 


23.7±0.2 


2002-2005 


z 


SDSS DR5 


21.2±0.3 


2002-09-05 




Co-add (DRSN1, DRsup) 


20.9±0.1 


2002-2005 


Z NTT 


NTT/EMMI 


20.69±0.06 


2006-07-18 


Y 


UKIDSS 


20.40±0.12 


2005-11-26 


J 


UKIDSS 


19.98±0.10 


2005-11-26 


K 


UKIDSS 


19.21±0.08 


2005-12-09 



4 CANDIDATE FOLLOW-UP OBSERVATIONS 

4.1 ESO NTT imaging observations 

Candidates were observed with filter [J611 in the ESO Multi- 
Mode Instrument (EMMI, Dekker, Delabre & Dodorico 
1986), at the 3.6-m New Technology Telescope (NTTfl, 
for 900 s on the night beginning 2006 August 18. The fil- 
ter/detector combination is referred to here as zntt and is 
similar to the SDSS z-band. Using the measured CCD sen- 
sitivity curve, and the filter transmission curve, we followed 
the procedure used by Hewett et al. (2006) to establish the 
colour relation zntt = z — 0.05(i — z) (in AB magnitudes), 
for dwarf stars, O to M. Photometry in i and z of SDSS DR5 
sources in the field was converted to zntt, and the bright- 
ness of the source was measured via relative photometry 
using a fixed aperture size. For a typical quasar at z = 5.9 
the difference between zntt and z is at most 0.05 mag. 

With the zjvTT-band photometry eight of the 10 ob- 
served objects were too blue in i — z and too red in z — J to 
be considered as likely quasars (see Table A.l for the names 
of these rejected candidates) . Such objects appear to be the 
main contaminant in our search. Due to photometric errors, 
the z-band magnitudes are scattered to a brighter value, 
making the objects appear redder than i — z > 2.2 and blue 
in z — J. 

The two remaining candidates, where the zntt magni- 
tude confirmed the co-added SDSS z-band magnitude, were 
selected for spectroscopic follow-up. 

4.2 ESO VLT spectroscopic observations 

Spectra of the two candidates were obtained on 2006 
September 1 and 2 using the FOcal Reducer/ low dispersion 
Spectrograph 2 (FORS2, Appenzeller et al. 1998) on the 
8.2-m Very Large Telescope (VLT) AntrQ. The candidates 
were observed through the 600z holographic grism with a 
l'.'O wide long slit. Individual exposures were 600 s dura- 
tion. The pixels were 2x2 binned to decrease the readout 
time and noise, giving a spatial scale Cf.'25 pixel -1 and a 
dispersion of 1.62 A pixel -1 . The nights were photometric 
with seeing around 0'.'5. For the wavelength calibration, ex- 
posures of He, HgCd and Ne arc lamps were obtained. The 

1 Based on observations made with ESO Telescopes at the La 
Silla Observatory under programme ID 077.A-0807 

2 Based on observations made with ESO Telescopes at the 
Paranal Observatory under programme ID 077.A-0310 
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Figure 2. Discovery spectrum of the z = 5.86 quasar ULAS 
J0203+0012. Below the object a noise and a sky spectrum are 
plotted, both in units of 10~ 17 ergs" 1 cm -2 A -1 . The spectra 
are boxcar averaged over three pixels. The two-dimensional spec- 
trum is shown at the bottom. The Lya (1216 A) and N V (1240 A) 
emission lines are marked. The object shows the characteristic 
continuum decrement across the Lya emission line, with very lit- 
tle flux blueward of the break. Around 8500 A a C IV AA1548,1551 
absorption doublet can be seen in the spectrum. The 'absorption 
lines' at 8917 A, 9315 A and 9518 A arc sky subtraction artifacts, 
and the absorption near 9073 A is caused by a bad pixel. 



rms of the wavelength calibration was better than 0.2 A. 
The spectra cover the wavelength range AA7370-10200. The 
spectrophotometric standard star GD50 (Oke 1990) was 
observed with a 5arcsec slit for relative flux calibration, 
and the absolute calibration was determined using the ob- 
served Y-band flux. Total exposure times were 1200 s for 
ULAS J020332. 38+001229.2 (hereafter ULAS J0203+0012) 
and 2400 s for ULAS J022905. 11+005506.5. 



5 RESULTS 

Fig. [2] shows the spectrum of ULAS J0203+0012. The spec- 
trum is characterized by a broad, asymmetric emission 
line and a continuum break around 8350 A. We interpret 
these features as the Lya emission line and the continuum 
break resulting from Lya forest absorption at a redshift of 
z ~ 5.86. The N V 1240 emission line is blended with the Lya 
line. The C IV 1549 emission line would lie at ~10630 A, out- 
side the wavelength range of the spectrum. We determined 
the redshift of the quasar by fitting a Gaussian to the red 
part of the Lya line. This gives a redshift of z =5.864+0.003. 
The peak of the emission is at 8342 A giving a redshift of 




Figure 3. Finding chart for ULAS J0203+0012. The chart is 
taken from the Y-band image and is 5'x5' in size. The image 
is boxcar averaged over 3x3 pixels. The quasar is marked by a 
cross-hair. 




7000 8000 10000 15000 20000 

Wavelength (A) 

Figure 4. Spectral energy distribution of ULAS J0203+0012 
(solid line) and a composite of three SDSS quasars with 5.8 < 
z < 6.0. 



2 = 5.862. The rest-frame equivalent width of the Lya + N V 
emission line is 64 A (uncorrected for Lya forest absorp- 
tion), which is comparable to other z ~ 6 quasars (e.g. Fan 
et al. 2004). A finding chart of the quasar is shown in Fig. 
[3] and the magnitudes are given in Table 1. It should be 
stressed that ULAS J0203+0012 is too faint in the optical 
passbands to be present in the SDSS DR5 catalogue. The 3<r 
2-band magnitude listed in Table 1 is derived using aperture 
photometry on the SDSS pixel image at the location of the 
UKIDSS source. 

The second candidate, ULAS J022905. 12+005506. 5, 
turned out to be a star, probably an M dwarf. 



6 DISCUSSION 

In Fig. [4] we compare the optical and infrared magnitudes 
of ULAS J0203+0012 with a composite of three 5.8 < z < 
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6.0 SDSS quasars which were observed during the science 
verification of WFCAM. Both spectral energy distributions 
(SEDs) are scaled to F v = 1 at the centre of the Y-band 
at 10305 A. ULAS J0203+0012 is significantly redder in the 
infrared than the SDSS quasar composite. If we approximate 
the SEDs in the infrared by a power law with slope a (F v oc 
v a ) then the slopes for ULAS J0203+0012 are a = -1.44 ± 
0.17 between the Y and K and a = —1.25 ± 0.21 between 
J and K. The SDSS quasars have average slopes of a = 
-0.52 ± 0.09 and a = -0.56 ± 0.12 between Y and K and 
J and K respectively. This is consistent with the average 
slope of a — —0.57 that was found for a sample of 45 SDSS 
quasars with 3.6 < z < 5.0 by Pentericci et al. (2003). 

ULAS J0203+0012 is significantly fainter than the 
SDSS quasars. The absolute magnitude at 1450 A in the 
rest-frame is M1450 = —26.2, nearly a magnitude fainter 
than that of the SDSS composite (M1450 = —27.1). How- 
ever, because ULAS J0203+0012 is redder, the extrapolated 
absolute £?-band magnitude is similar to that of the SDSS 
composite, Mb = —27.9 and Mb = —27.7 respectively. 



7 CONCLUSIONS 

We have reported initial results of a survey program us- 
ing UKIDSS LAS data to discover high-redshift quasars. 
Our results from 106 deg 2 demonstrate that the federated 
UKIDSS and SDSS data can be used to carry out such a sur- 
vey. The analysis and follow-up strategy we have described 
is relatively efficient and should be applicable to the entire 
UKIDSS LAS survey area of 4000 deg 2 . 

Follow-up observations of the UKIDSS DR1 quasar 
candidates in the SDSS Stripe 82 region are now com- 
plete. It is difficult to draw any strong conclusions regard- 
ing the space density of quasars from one object but we 
can estimate the number of quasars we would expect in 
the UKIDSS survey by extrapolating the z — 6 luminos- 
ity function of Fan et al. (2004) to fainter magnitudes. The 
density of high-redshift quasars as found in the SDSS is 
<E>(Mi450 < -26.7) = 6 x 10~ 10 Mpc~ 3 (Fan et al. 2004) at a 
redshift of z ~ 6.1. We assume that the comoving density of 
quasars evolves with redshift as log($)oc — 0.43z (Schmidt, 
Schneider & Gunn 1995). 

Number density predictions are calculated for the Y- 
band, corresponding to a rest- frame wavelength of ~ 1450 A 
at z — 6.1, for which no correction for the presence of Lya 
emission or the Lya forest is required. Integrating down to 
Yab = 20.4 (the median 7a limit of the DR1 data) we expect 
one quasar in ~130 deg 2 in the redshift range 5.8 < z < 6.5. 
This is consistent with the one object, ULAS J0203+0012, 
that was found in 106 deg 2 . We can estimate in the same 
way the surface density of z > 6.5 quasars in the LAS and 
we find one object per ~390 deg 2 down to Yab = 20.4. The 
UKIDSS LAS two-year plan (Dye et al. 2006) will cover 
around 2100 deg 2 . We expect this area to contain ~ 22 new 
z > 5.8 quasars (excluding six quasars already discovered 
by SDSS), of which about six would be at z > 6.5. 

Above a redshift of z ~ 7.5 the Lya line shifts out of the 
Y-band. Potentially quasars above this redshift could be dis- 
covered as Y-band dropout objects. However, the predicted 
number density for such quasars is very low (< 1 quasar in 



Table Al. Quasar candidates that were rejected after the follow- 
up observations. Using the zjvTT _ band photometry (Sect. 4.1) 
these eight objects are too blue in i — z (i — z < 2.2) and too red 
in 2 — J to be considered as likely quasars. 



UKIDSS name 



ULAS J030352. 
ULAS J003852. 
ULAS J022119. 
ULAS J223845. 
ULAS J021753. 
ULAS J223457. 
ULAS J011650. 
ULAS J021047. 



55+003334.1 
24-005324.0 
09-010905.9 
14+011132.8 
93+003906.9 
11+005136.6 
20+011532.2 
00-001419.6 
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